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. Abstract. I review current theories of disk galaxy formation. All fail in 

the light of recent data for low surface brightness galaxies. The difficulty 
stems from the presence of a unique acceleration scale in the data. 
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All we know so far is what doesn't work. 

— Richard Feynman 
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C"*** . 1. Introduction 

<N 

Galaxy formation is an intriguing but difficult subject. Substantial progress 

■ has been made over the last few years, culminating in a rough consensus on a 
OO ! "standard" picture of galaxy formation (e.g., Mo et al. 1998). In spite of my 

■ own modest contributions to this (e.g., Mo et al. 1994), it is my contention that 
^ cj this picture is wrong. 

Qh' Empirically, there are a number of striking regularities in the dynamical 

q . data for disk galaxies. These enable us to describe the effective force law. It is 

5h ' remarkably uniform, and depends only on the luminous mass. 



2. Galaxy Formation 



^ ' The general picture we have of disk galaxy formation is of cool gas settling into a 

disk in a potential well dominated by dark matter. The dark matter is assumed 
to reside in a dissipationless, effectively spherical distribution, and to interact 
with the potentially luminous ( "baryonic" ) matter only gravitationally. To make 
predictions from this general picture, we need to specify something about both 
the halos and the baryons. I would like to emphasize the importance of making 
a priori predictions. If we simply tune our favorite model to fit to the data, we 
will certainly fit the data. 

At minimum, two parameters are required to describe halos, such as a 
characteristic mass Aih and size Rh- Further parameters describing axial and 
kinematic isotropy, etc., could be invoked, but these distract from the basic point 
and are not fundamental. For the baryons, many parameters might in principle 
come into play: initial angular momentum and gas temperature, conversion of 
gas into stars and feedback into the ISM, and so on. It would be nice to study 
a population of galaxies for which the effects of the baryons were minimized. 
Low surface brightness (LSB) galaxies, with central surface brightnesses (jlq > 
23 B mag. arcsec" 2 , turn out to be just such objects. They occupy a different 
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region of parameter space than do high surface brightness (HSB) spirals, so they 
provide genuinely new tests of ideas contrived to explain observations of HSB 
disks. 

2.1. Density Begets Density 

Since I am going to contradict some widely held notions, it seems only fair to 
start by trashing my own favorite idea. Many of the properties of LSB galaxies 
(McGaugh 1992) suggested to me that they were basically stretched out, lower 
density versions of HSB galaxies. They also appear to be somewhat younger than 
HSB galaxies, suggesting a later collapse epoch (McGaugh & Bothun 1994). 

It seems natural that the properties of dark halos might dictate the prop- 
erties of the luminous galaxies they contain. In addition to a distribution of 
masses Mh which gives rise to the luminosity function, there could also be a 
distribution of scale sizes Rh which, at a given mass, gives rise to the distribu- 
tion of disk scale lengths (or equivalently, the surface brightness distribution: 
McGaugh 1996). The mapping from (Mh, Rh) to (L, h) need not be simple, but 
assuming it is forms an obvious starting point. In this picture (which I refer to as 
'density begets density' or DD), HSB disks arise from large density fluctuations 
<5, and LSB disks from low 5 (but not necessarily low mass). 

DD makes several predictions beyond the data which motivated it (Mc- 
Gaugh & de Blok 1998a). One is that since HSB and LSB galaxies represent 
populations which arose from different characteristic 5, there should be a shift 
between the correlation function of the two. This prediction was confirmed (Mo 
et al. 1994). 

In DD, LSB galaxies and their halos are stretched out versions of HSB 
systems. Two predictions follow about their dynamics. They should have slowly 
rising rotation curves (which is true), and they should deviate systematically 
from the Tully-Fisher (TF) relation. 

The expected departure from the TF relation goes in the sense that LSB 
galaxies should rotate slowly for their luminosity. This follows simply because 
Rh is large for a given A4h, and V 2 = GAi/R. There is some freedom to tune 
the amount of the shift, so long as it is systematic with surface brightness. It 
would not fit nicely into this picture if there were zero shift. 

This is exactly what is observed (Figure 1). Regardless of surface brightness, 
LSB galaxies fall on the sameQ TF relation with the same normalization as HSB 
galaxies. The shift I expected is illustrated by the lines in Figure 1, and clearly 
does not occur. 

2.2. Same Halo 

As an alternative to DD, I have investigated the class of theories which I gener- 
ically label 'same halo' (SH; McGaugh & de Blok 1998a). The basic idea here is 

1 This result has been independently obtained by Zwaan et al. (1995), Sprayberry et al. (1995), 
Hoffman et al. (1996), and Tully & Verheijen (1997). It has been disputed by Matthews et al. 
(1998), who suggest that some faint LSB galaxies fall systematically below the TF relation. 
This is a result of comparing to a TF relation with a shallow slope determined from a fit to 
bright galaxies. The data of Matthews et al. are consistent with other data (Fig. 1). Indeed, 
they add weight to the faint end which is much needed for constraining the slope. 
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Figure 1. The luminosity-rotation velocity (Tully-Fisher) relation. 
Data are shown for galaxies with V c measured from the flat part of 
resolved, extended rotation curves (solid symbols; de Blok et al. 1996) 
and for some with recent high quality velocity widths (open symbols; 
Matthews et al. 1998). Symbols distinguish different bins of central 
surface brightness — stars: /x < 22; squares: 22 < /i < 23; triangles: 
23 < no < 24; filled circles: fio > 24. Galaxies with open symbols 
predominantly have > 23. The data do not distinguish themselves 
by surface brightness in this diagram. The lines are not fits to the data, 
but have the "virial" slope of —10. They are labeled by a correspond- 
ing central surface brightness illustrating the systematic shift that is 
expected from simple arguments (see text). 



that there is no distribution in R^. Simulations of halo formation (e.g., Navarro 
et al. 1997; hereafter NFW) indicate a strong correlation between halo parame- 
ters, consistent with this picture. At a given mass, the halo is the same regardless 
of the scale length of the optical galaxy. This yields the desired TF relation, by 
construction. 

Since there is now no distribution in R^, we must invoke some other mech- 
anism to give the observed distribution of optical surface brightnesses. This 
is usually assumed to follow from the initial angular momentum. In terms of 
Peebles's spin parameter A, the scale length of the luminous disk is h ~ \Rh- 
The precise equation can be more complicated (Dalcanton et al. 1997Q Mo et 
al. 1998), but this encapsulates the basic idea. Baryons in a halo with low initial 



Given the emphasis placed on a priori predicti ons, it is worth no ting that the original version 
of the scenario described by Dalcanton et al. ( istro-ph/9503093) predicted a shift in the TF 
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spin collapse a long way before rotational support is achieved, forming an HSB 
galaxy with a short scale length. A high spin halo of the same mass forms an 
LSB galaxy with a much larger disk scale length. This makes the rather dubious 
assumption that there is no interaction between disk and halo which can transfer 
angular momentum between the two. 

The surface brightness distribution is now determined by the initial distri- 
bution of A rather than 6, causing a different problem to arise. In fixing the 
failings of DD with regards to the TF relation, we lose its success in predicting 
the shift in the correlation function. Simulations show no correlation between 
spin and environment (e.g., Barnes & Efstathiou 1987). Therefore there should 
be no shift in the correlation function with surface brightness as is observed. 
This is as much a problem for the SH picture as the TF relation is for DD. We 
might or might not be able to fix it (cf. Mihos et al. 1997; Moore et al. 1999), 
but whatever we come up with is a patch after and against the original fact. 

2.3. Mass and Light 

There is information in the rotation curve beyond the TF relation. The shapes 
of rotation curves are also related to luminosity, as noted by Rubin et al. (1985) 
and Persic & Salucci (1991). Though I would not claim as strict a relation as 
implied by the 'universal rotation curve' of Persic & Salucci, a correlation does 
exist and provides an additional test (Figure 2). 

LSB galaxies adhere to the relation between luminosity and rotation curve 
shape, provided that the radius is measured in units of the disk scale length 
(see also Verheijen, these proceedings). This implies that a good estimate of 
the rotation curve of any galaxy can be made from measurements of only two 
photometric parameters (L,h). Even though the dynamics are dominated by 
dark matter, we need only know the distribution of luminous matter to predict 
the rotation curve. 

This strong coupling of mass and light (as long stressed by Sancisi) is a 
general problem. It always leads to fine-tuning paradoxes, with the tail wagging 
the dog. The dominant, spherical halo composed of non-baryonic dark matter 
simply should not be so intimately related to the details of the luminous disk. 

2.4. CDM and NFW 

I have stressed the importance of a priori predictions. There has been sub- 
stantial progress in numerical investigations of the formation of halos in CDM 
simulations (e.g., NFW). These result in an apparently universal profile describ- 
able by two parameters, a concentration c and a scale V2oo- These are tightly 
correlated, effectively forming a one parameter family for any given cosmology. 

The data for LSB galaxies provide a good test of the these predictions 
(Figure 3). Because they have such large mass discrepancies, complications due 
to the baryonic component are minimized. The precise value of Y*, adiabatic 
compression, feedback, etc., simply do not matter. 



relation with surface brightness. The observed lack of such a shift was reported at about this 
time. Zero shift is 'predicted' in the published version. 
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Figure 2. The systematic dependence of the shape of rotation curves 
on luminosity and disk scale length. R34 is the radius at which V{R) 
has reached 3/4 of the asymptotic flat velocity. A large -R34 corresponds 
to a slowly rising rotation curve. R34 varies greatly between galaxies 
of the same luminosity but different surface brightness when R34 is 
measured in kpc. However, there is a good correlation in this diagram 
when i?34 is normalized by the disk scale length h. Some, and perhaps 
most, of the scatter is attributable to observational uncertainty. Also 
shown are model predictions of Dalcanton et al. (1997), with numbers 
labeling points by the logarithm of the halo mass. The models track in 
the opposite sense of the data, a problem generic to SH models which 
have not been tuned to fit the data. 



The predicted shape of the inner profile of the dark matter halo goes as 
p oc r -7 with 7 = 1 according to NFW. Moreover, the halo parameters are 
specified by the cosmology. Figure 3 shows, in several cosmologies, the expected 
rotation curves for NFW halos chosen to be a close match to the illustrated 
galaxy. The predictions fail in two ways: NFW halos have rotation curves 
which rise too steeply, and require very low disk-to-total mass ratios. 

The first failing is obvious by inspection, as first noted by Flores & Primack 
(1994) and Moore (1994). The data for LSB galaxies confirm and extend the 
conclusions of these works. There have been various attempts to wiggle out of 
this problem, without success (McGaugh & de Blok 1998a). Adiabatic contrac- 
tion should occur at some level, but acts in the wrong direction. One might 
suppose instead that dwarf and LSB galaxies suffered massive baryonic blow 
outs following an episode of intense star formation which rearrange the dark 
matter in the required way (Navarro et al. 1996). There is no empirical evidence 
that this ever happened. Gas has not been blown away (McGaugh & de Blok 
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Figure 3. The rotation curve of the LSB galaxy F583-1. The solid 
points are the data; the open points are the remainder after the sub- 
traction of the maximum disk component. This illustrates the large 
mass discrepancy typical of LSB galaxies, even at small radii. Also 
shown are the rotation curves due to NFW halos expected in several 
cosmologies. These rise too steeply at small radii, and require very low 
disk-to-total masses (/&) in order to match the velocity at large radii. 



1997), nor is there any debris around LSB galaxies (Bothun et al. 1994). It now 
appears that the idea is unworkable even in theory (MacLow &; Ferrara 1998). 

The problem is soluble if the amount of dark matter in the inner parts of 
the initial halos can be made low enough. In the context of NFW halos, this 
means very low c. This in turn requires a very contrived cosmology (Figure 
4). Standard (f2 = 1) CDM produces halos which are much too concentrated. 
Low f2 models (with or without a cosmological constant) fare better, but still 
produce halos in which c significantly exceeds the upper limits imposed by the 
LSB galaxy data. Using Navarros's code, the only way I have found to lower c 
further is to either reduce the normalization of the power spectrum well below 
the values inferred from COBE or clusters of galaxies, or to choose cosmological 
parameters (fi, Hq) which yield a shape parameter T < 0.1. This is inconsistent 
with the large scale structure constraint 0.2 < T < 0.3 (Peacock & Dodds 1994). 

This is very bad for CDM, or at least the NFW realization thereof. Could 
something be wrong with the profile suggested by these simulations? Initially, 
there seemed to be a good consensus between different simulations. More re- 
cently, this has become rather controversial. Most workers give 7 = 1 or 7 > 1 
(Dubinski 1994; Cole & Lacey 1996; Moore et al. 1998; Tissera & Dominguez- 
Tenreiro 1998; Nusser & Sheth 1998). Any 7 > 1 is inconsistent with the data. 
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Figure 4. Predicted and measured NFW halo parameters (c, V-joo)- 
These are expected to be tightly correlated in any given cosmology 
(lines). Different cosmologies are represented by different lines, as la- 
beled. The data place upper limits on c by assuming stars have no 
mass (Pickering et al. 1997). This is unrealistic in HSB galaxies, lead- 
ing to uninteresting limits (dots). The stellar mass is negligible in LSB 
galaxies, which place strong upper limits (crosses with arrows) on c. 
No plausible cosmology produces sufficiently low concentrations c to 
satisfy the data, even with the favorable assumption of low H$. 

On the other hand, Kratzsov et al. (1998) give 7 ~ 0.2. This is inconsis- 
tent with other simulations but consistent with the data (see Primack, these 
proceedings). Indeed, Kratzsov et al. (1998) constrain 7 by fitting the data for 
dwarf and LSB galaxies. This constrained profile may be consistent with their 
simulations, but it hardly constitutes an independent test. 

There is another problem in the baryon fraction. The general picture is 
of a universal baryon fraction /{, which is constant from halo to halo (giving 
rise, among other things, to the TF relation). A great deal of emphasis has 
been placed on measurements of in clusters of galaxies (e.g., Evrard 1997). 
However, the same exercise can be done in other systems as well, with the result 
that there is no indication of a universal baryon fraction (Figure 5). Indeed, 
a number of limits are in blatant contradiction. There are various outs, too 
numerous to discuss here (see McGaugh &: de Blok 1998a). To reconcile the 
cluster value of 0.1 with the LSB galaxy value of fb ~ 0.02 suggested in 
Figure 3, one might suggest that not all of the baryons have been incorporated 
into the disk. If so, there is no reason to think the fraction which have been will 
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Figure 5. Measurements of and limits on the baryon fraction in var- 
ious systems (McGaugh & de Blok 1998a). A single universal baryon 
fraction is not evident. Indeed, a number of limits are in conflict. For 
example, the limits for LSB galaxies refer only to the edge of the ob- 
served rotation curve (de Blok & McGaugh 1997). If we analyze them 
in the same way as clusters, the limits at lower left become points in 
the upper left portion of the diagram. 



be the same for all galaxies. Yet it must be. Any scatter would propagate into 
the TF relation, for which there is no room in the error budget. 

I have long been a believer in the CDM paradigm. Yet I am forced to 
conclude that the accumulated evidence now weighs heavily against it. Unfor- 
tunately, the existence of dark matter is not explicitly falsifiable. Theories in 
which the dominant mass component is invisible have sufficient flexibility to 
accommodate anything. 



3. A Physical Scale in the Mass Discrepancy 

Galaxy formation theory fails to predict a mass distribution which yields the 
correct effective force law. So let us take a more empirical approach. A useful 
measure of the severity of the mass discrepancy is the ratio of the required 
gravitating mass to the total observed luminous mass. One would not expect 
dark matter to be aware of any particular physical scale. Is there one in the 
data? 

For some reason, our brains think first in terms of linear size. Galaxies 
are big. Is there anything special about the mass discrepancy at some large 
length scale? The answer is a resounding no (Figure 6a). There are galaxies for 
which the mass discrepancy is not apparent until quite large radii, and others 
(predominantly LSB galaxies) in which it appears nearly at R = 0. Size does 
not matter. 
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Figure 6. The mass discrepancy as a function of (a) radius and (b) 
orbital frequency (uj = V/R). Each point represents one resolved 
measurement in the rotation curve of a disk galaxy. Data for many 
galaxies of all luminosities and surface brightnesses are plotted to- 
gether. The mass discrepancy is defined as the ratio of dynamical 



mass Mtot = RV /G to luminous mass M 



lum 



stellar mass-to-light ratio Y* = 2 is assumed, as is M 



T*L + M gas . A 
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A line of unity indicates no mass discrepancy. There is no preferred 
scale in either plot. The mass discrepancy can appear at small as well 
as large radii (a), where the eye can perceive some individual rotation 
curves. This is less true in (b), which shows that the mass discrepancy 
occurs preferentially at low frequencies, though there is a lot of scatter 
at a given uj. Data are taken from the compilation of Sanders (1996) 
and de Blok & McGaugh (1998). 



There are other scales besides linear size. A plot of the mass discrepancy 
against orbital frequency (Figure 6b) begins to show some organization. While 
there is still a lot of scatter, this is a hint that there may be some interesting 
scale. 

A scale which is truly unique to galaxies is that of low acceleration. The 
centripetal acceleration which keeps stars in their orbits is typically < 1 As -2 , 
only one part in 10 11 of what we experience at the surface of the earth. The 
data show an enormous regularity in terms of acceleration (Figure 7: see also 
Sanders 1990). It matters little how we compute the luminous mass (assuming 
constant Y* or maximum disk). The reason for this is that LSB galaxies have 
large mass discrepancies in either case, and also have the lowest accelerations. 

To the best of my knowledge, there is no dark matter based prediction of 
this phenomenology in the literature. Attempts to impose it by fixing the surface 
density of dark halos fail (Sanders & Begeman 1994). It makes no sense to me 
that dark matter should be aware of a particular physical scale in this fashion. 
This is reflected in the difficulties discussed in §2. 

The regularity of the data strongly suggest that some universal phenomenon 
is at work. There is one suggestion that the usual dynamics should be modified 
at an acceleration scale: MOND (Milgrom 1983a). Milgrom (1983b) made a 
series of predictions specific to LSB galaxies. I find this quite remarkable, as at 



9 





a (m s ) a (m s ) 



Figure 7. Like the preceding Figure, but now with the abscissa in 
terms of centripetal acceleration V 2 /R. The data all fall together: it is 
no longer possible to distinguish individual rotation curves. The mass 
discrepancy in spiral galaxies occurs uniformly at a particular acceler- 
ation scale of ~ 10~ 10 ms~ 2 . This occurs regardless of what we assume 
for the mass-to-light ratio, as the mass discrepancy is large in LSB 
galaxies even in the maximum disk case (b). The only way to avoid 
the appearance of a particular acceleration scale is to allow T* — > 0, 
in which case there is a mass discrepancy everywhere. Moreover, it 
is possible to choose an optimal mass-to-light ratio which minimizes 
the scatter (c) as far as possible given experimental uncertainty. The 
value of each galaxy's Y* in (c) could be chosen to minimize the scatter 
in this plot, but this is not what has been done. Instead, a specific, 
universal prescription is applied: the optimal value of Y* is that indi- 
cated by MOND fits to the rotation curves. MOND resolves the mass 
discrepancy in disk galaxies (d), providing an accurate description of 
the effective force law therein. 
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that time much of the debate over LSB galaxies was about whether they existed 
at all, with the majority opinion being negative. Now we know they exist, and 
have the data to test Milgrom's predictions. 

Every prediction Milgrom made about LSB galaxies in 1983 is confirmed 
(McGaugh & de Blok 1998b). 

There are many claims to have falsified MOND. An extensive survey of the 
literature turns up no credible empirical evidence against Milgrom's hypothesis 
(McGaugh & de Blok 1998b). Indeed, there are many cases which are cited as 
evidence against MOND which can equally well be argued to support it (dwarf 
Spheroidals being one example). Are we failing to see the forest for the trees? 

It remains an open question whether MOND works in all places where a 
mass discrepancy is inferred. But it is extremely effective at resolving the mass 
discrepancy in disk galaxies (Figure 7d). This is telling us something. For all 
practical purposes, MOND is the effective force law in disk galaxies. 

Irrespective of whether MOND is correct as a theory, it does constitute an 
observed phenomenology. As such, it provides a very strong test of dark matter 
galaxy formation theories. The requirement is this: for any model, it must be 
possible to apply the MOND formula to the disk component and obtain the 
rotation curve which is actually produced by the combination of disk and halo 
components. 

There is very little freedom to achieve this. All the many parameters of 
conventional models must be encapsulated in one tightly constrained parameter, 
Q. This is the ratio of the stellar mass-to-light ratio that is required to obtain 
a MOND fit to that which is 'correct' in the model (Q = Tf OArD /T? M ). g ince 
the mass-to-light ratios indicated by MOND are already quite reasonable for 
stellar populations (Sanders 1996, de Blok & McGaugh 1998), this is a very 
restrictive requirement indeed: Q ~ 1 with no other adjustable parameters. 
Nevertheless, this is what is required to reproduce the observed phenomenology. 

Until there is a credible explanation for the MOND phenomenology in the 
framework of the standard paradigm, we should be as skeptical of the existence 
of dark matter as we are of the need to modify Newton's Laws. 

Acknowledgments. I am grateful to Erwin de Blok for our long and fruit- 
ful collaboration, and to Vera Rubin, Renzo Sancisi, Bob Sanders, and Jerry 
Sellwood for many animated discussions of disk dynamics. 

Discussion 

Spergei. I assert that an acceleration scale does occur naturally in dark matter. 
McGaugh: Did you predict this in 1983, or did Milgrom? 
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